Hysteretic Depinning of Anisotropic CDW 
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We investigate the depinning transition in a dirty periodic medium considering a model of layered 
charge density waves as a prototype system. We find that depinning from strong disorder occurs 
via a two stage process, where, first, the pinned system experiences a continuous transition into a 
plastically sliding state and undergoes a second sharp hysteretic transition into a coherently moving 
3D state at higher drives. In the weakly disordered system the depinning into a coherently sliding 
state remains continuous. 



Depinning of periodic structures such as charge den- 
sity waves (CDW), vortex and domain wall lattices, and 
Wigner crystal from a random pinning potential under 
the influence of an external driving force is one of the 
paradigms of condensed matter physics. All these sys- 
tems share one thing in common: many of elastically cou- 
pled degrees of freedom interacting with a quenched ran- 
dom environment. In this Letter we report analytic re- 
sults on the nature of the depinning transition in dirty pe- 
riodic media using CDW dynamics as a prototype model 
and discuss possible extensions to other systems. 

Two types of depinning have been observed: a smooth 
non-hysteretic transition with a unique pinning thresh- 
old, and transport switching characterized by an abrupt 
hysteretic transition into a sliding state H, Smooth 
depinning described in terms of critical behavior || fol- 
lows from the description of the above systems as clas- 
sical field associated with the distortions of the system. 
Much of the switching behavior is explained by the pos- 
sibility of plastic deformations allowing the amplitude of 
CDW to vanish along certain surfaces within the system 
||,[|]. Recent experimental and numerical studies |^-||] 
of vortex transport in HTS have demonstrated that the 
depinning of the vortex lattice can also be accompanied 
by plastic effects. 

These latter findings suggest that plastic effects can 
play an essential role in the depinning, and that the 
nonequilibrium steady state near the transition resembles 
fluid-like motion. On the contrary, at very high veloci- 
ties well above the depinning transition the influence of 
disorder on the dynamics is suppressed and one can ex- 
pect coherent motion of an almost perfect solid periodic 
system. A problem of separation of these two different 
driven regimes has been addressed in |J in the context 
of vortex transport. It was proposed that the driven pe- 
riodic medium subject to sufficiently strong disorder un- 
dergoes a sharp hysteretic dynamic transition from co- 
herent motion with almost perfect structure to fluid-like 
plastic dynamics upon decreasing the drive at a second 
critical force Ff force well above the pinning threshold 
Ft. This was called dynamic freezing. Upon increasing 



the driving force from the pinned state the vortex lattice 
starts to slide at F = Ft, this depinning being followed 
by the multiple plastic effects, and elastic motion recovers 
at F — Ff > Ft- This concept received strong support 
both from earlier observations of plastic effects Q and 
from the subsequent transport measurements on MoGe 
superconducting films fiof] . 

The prediction of the possibility of a dynamical phase 
transition in the driven state B was later expanded by 
nice scaling arguments onto three dimensional CDWs 
|n| . The properties of the driven coherent phase were 
examined and discussed in [p~2| — p~6|| - Yet a number of 
unresolved problems remains. The fundamental issue is 
the nature of the depinning transition (continuous vs. 
switching), and the question under what conditions ei- 
ther type occurs. How would dynamic freezing evolve 
with decreasing strength of the disorder? What are the 
conditions for the existence of the plastic flow regime as 
the state intermediate between the pinned and coherently 
moving states? Is it possible to have depinning directly 
to a coherently moving state? 

In this Letter we address these questions using CDW 
transport as an example system. We consider a model of 
an anisotropic CDW [|l9| and develop a nonperturbative 
self-consistent description of the dynamic transitions in 
driven dirty systems. We find that if disorder is suffi- 
ciently strong, depinning occurs in two stages: first the 
CDW depins in a driven decoupled state where 2D CDWs 
in each layer slide independently, and second, upon fur- 
ther increasing the driving force the system experiences 
a second transition into a coupled coherently moving 3D 
phase. This "sequential depinning," corresponds to a dy- 
namic freezing transition scenario proposed in |^| , where 
the periodic system first depins into a plastically moving 
state and then upon further increase of the drive expe- 
riences a transition into a coherently moving dynamic 
state. In the system with weak disorder decoupling does 
not occur. 

The paper is organized as follows. First we describe 
the model and derive a selfconsistent equation for the 
shear modulus. Then we analyze the cases of weak and 
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strong pinning, and determine the disorder induced dy- 
namic decoupling (the instability point) as the current 
at which the onset / disappearance of the coupling occurs. 
In conclusion we construct the general dynamic phase 
diagram for periodic structures driven through quenched 
disorder. 

The overdamped dynamics of an anisotropic layered 
CDW are governed by the equation (CDW moves along 
the z-direction): 

Xfa(x, t) = 7V 2 0i(x, t) + 7/i [sin(^j+i - fa)+ 

sin(0i_i - fa)] + F + ~fVsixi(fa - a<(x)). (1) 

Here <fi is a CDW phase (displacement field), A is a fric- 
tion coefficient, v is the average velocity of CDW, 7 is the 
elastic constant, ho is the anisotropy parameter charac- 
terizing layer coupling, V is the strength of the random 
potential, a(x) is a random phase, i is a layer index, and 
x is a D-dimensional vector in the layer. In the CDW 
models a comes from the backscattering part of the im- 
purity potential: a;(x) = 2ki?(x + Zo«a), where kp is the 
Fermi vector, and r = (x, zoia) is the coordinate of the 
impurity (zo is the unit vector in the z direction). Since 
the in-plane pinning correlation length exceeds the impu- 
rity spacing the random phase a is commonly considered 
as a random variable homogeneously distributed in the 
interval [0, 2tt] @. 

If the anisotropy parameter ho is large one recovers 
the continuous limit, then i — » z and the term in square 
parenthesis becomes simply d 2 (f>/dz 2 . The related Hamil- 
tonian has the form: 



n = 7 / d° x — 

a 



~ (Vfa) 2 +~Ho(d z <P) 2 +V(x,z, 



(2) 



force F Ctd = 7 (F 4 /Mo) 
correlation length is £ c ,d = (V 4 /ho) d , d = D 



where we set a = 1. This D + 1-dimcnsional system 
experiences a continuous depinning transition at critical 

l/(4~d) -, , . . 

Ho ^S> 1, and the pinning 
1. 

If /io — ► the decoupling transition occurs in a static 3D 
system at ho < Hmin — V 2 . Note that in a decoupled 
limit the equation of motion ([!]) becomes isotropic and 
Ho — > 1, d — ► D. 

To obtain a coarse-grained description in terms of a 
slowly varying < part of the phase we integrate out 
its fast component <j> = j- f Q ° tp(t)dt, to = l/v, I is the 
CDW period (hereafter we will drop the bar). The coarse 
graining procedure is straightforward in case of weak (as 
compared to coupling) disorder, ho 3> vV, where phase 
variations from layer to layer are small and sin(</>j_|_i — 
fa) & {4>i + \ — fa). Going over to a continuous description 
one arrives at the coarse-grained equation of motion for 
the slowly varying part of the phase: 

Xfax, z) = 7 (V 2 + nd 2 z ) - \vd x cb + F p (x, z), (3) 

where h ~ Ho an d Fp is the effective random force (ran- 
dom mobility-like term) originating from the random sin 



field (analogous to ]T3]) with the correlator defined as 
CFJx,z)F p (x',z')) = A(v)S(x - x')6(z - z') with g§, 

Hi 



A(v) = 




v > 7V/ I / A 



(4) 



We omitted the coarse-grained pinning potential term 
f p (4>) and the disorder-induced KPZ term: it can be 
shown that in the perturbative high velocity limit 
those terms are irrelevant. 

We are interested however in the strong disorder 
regime where the static CDW is decoupled. In this 
case the coarse-graining procedure is more subtle. To 
carry it out note that in the high-velocity case, v 3> 
7y / /I/A , A(v) does not depend on the elastic constants 
7 and the anisotropy parameter /i (provided the strength 
of disorder jV is fixed). This reflects the fact that at 
high velocities A(v) is controlled by the dissipative part 
1/iXv of the response function and enables us to carry 
out the coarse graining procedure also in the case of 
weak anisotropy Ho "C VV (i.e. when a replacement of 
sin(0j + i — by fa+\ ~~ 4>i is not possible). The coarse- 
grained equation of motion reads (v 3> 7 V / /I/A): 

Xfa(x) = 7V 2 (^(x) + 7/x [sin(0 i+ i - fa) 

+ sin(^_i - fa)] - \vd x 4>i + F p> i(x). (3') 

It is convenient to rewrite the Eq. (3') in the form 

^ G^fa = 7Mo [sin (fa+i - <f>i) + sin (fa-i - fa)] 
3 

-7/i [fa+i + fa-! - 2 fa] + F p>i (x) + £j(x) (5) 



with 



G^(t) = (\d t + Xvd x 



7 v 2 + 27/x)<y« 

-J/J>(Si,i+i + 



where £j(x) denotes a source term, which wil be sent 
to zero at the end of the calculations. One can now 
solve eq. (^) iteratively, generating an infinite number 
of tree diagrams with either F p ^ (x) or Ei (x) at the ends 
of each branch, and averaging subsequently over the ran- 
dom force keeping only the linear in £j(x) terms. To find 
the self-consistent equation for h we note that the self- 
consistency condition requires that the phase-containing 
terms in the r.h.s. of eq. (||) canceled each other. After 
that we obtain in the lowest order in ho'- 



H = Ho exp 



:((0(x, z,t)-fax, z + a,t)f) 



(6) 



Note also that since the correlations of F p ,i (x) are Gaus- 
sian the higher order cumulants do not appear within 
this scheme. 
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The correlation function in the exponent is calculated 
with the Hamiltonian from (||) with a = 1: 



C(v, fi, A{v)) = ((0(x, z, t) - </>(x, z + a, t)) 1 ) 
1 r d 3 k A(u)(l - cosfc z a) 

7 2 7 (27r) 3 ( fc 2 +k 2 + ^ fc 2) 2 + A2 W 2 fc 2/ 7 2 



(7) 



The main contribution comes from the maximal k z , 
therefore we can replace k z — > 7r/a in the integrand. 
Then in the 3D case one arrives at 



C(£, M ,A(«)) 



A(«) 



(4 7 r) 3 7 2 / « 3 / 2 



(8) 



where C = Aw/7. Outside the critical region u w -F/7, 
and within the critical region v ~ (F — J F c ) 1 ~( 4_d )/ 6 /7 
pfj . Note that in the limit of very large velocities Eq. (|7]) 
provides the expected behavior [i — > and the system 
always remains coupled. 

Now comes the central point of our discussion: solving 
the self-consistent equation for fi. The disappearance of 
the solution to Eq. (|jj) implies decoupling of the sys- 
tem. To capture the transition we will be seeking for the 
moment of the first disappearance of the solution. For 
the sake of simplicity we can replace the expression in 
square brackets in (7) by 27r 2 /x/ \/v 2 + 4-7r 2 /x, which gives 
the same asymptotics in the limit of small and large ve- 
locities. As we will shortly see in the case of strong dis- 
order the decoupling transition occurs at large velocities, 
whereas in the weak disorder case there is no decoupling. 

Strong disorder: V 2 > fi . In the limit of large veloci- 
ties Eq.(o) assumes the form: 



fj, = exp 



V 4 



32tt£V /2 



(9) 



To find the point /i c of the disappearance of the solu- 
tion we derivate both sides of (^|), obtaining the condition 
V 4 = 647to 3 /z 1 / 2 and find 



Mo/e 2 



V 4 



((> 1 \ //.- 



1/3 



(10) 



The last task to complete our calculation is to verify 
that the critical decoupling velocity indeed falls into a 
large velocities interval. To this end note that the large 
velocity condition v c 3> reduces at the instability 
point to V 2 ^> i^/TifjLo/e 2 which is just the condition of 
the strong disorder assumed and therefore our assump- 
tion is justified (one has to bear in mind that in our 
dimensionless units /1 <C V 2 <C 1). The critical velocity 
( |l0| ) agrees with the result of |Q suggested by nice scal- 
ing arguments if one substitutes A(w) from (||) instead of 
the unspecified mean squared pinning strength g of p| . 

To understand the meaning of the obtained result let 
us construct the F — v dependence starting with the as- 
cending branch. Below the 2D critical force F 2D the de- 
coupled system remains pinned. At F = F 2D the system 



undergoes a smooth depinning transition into the plasti- 
cally moving state in which the system remains decoupled 
and each layer moves independently. Upon further in- 
crease of the drive, the mean velocity of the CDW reaches 
the critical value and the system gets coupled into a 3D 
moving state. Note that since in the 3D regime pinning 
force experienced by the moving CDW, F 3D (v), is less 
then the corresponding pinning force in the 2D regime: 
F 3D (v) ~ F 2D {v)/^/Ll(v) < F 2D {v), where L±(v) is 
the (velocity dependent) correlation length across the 
layers, the pinning correction to the 3D velocity is smaller 
than the corresponding correction to the 2D velocity. As 
a result the 3D branch of the F — v dependence lies above 
the i>2d{F) curve and the transition from the plastic to 
the coupled elastic motion at v = v c upon increasing 
drive acquires an abrupt switching character (see Fig.l). 
Going down from the high velocities the system follows 
first the elastic 3D behavior with fi ps /xo, and then as 
velocity decreases to v — v c , the system decouples and 
jumps down to the 2D branch corresponding to the plas- 
tic motion (see Fig. la) at F = F c ^ own ) < F c r U p) ■ There- 
fore in the limit of strong pinning the transition from 
plastic to elastic motion is a switching hysteretic transi- 
tion. 

An immediate reservation regarding the proposed sce- 
nario is in order. The above picture suggests that the 
real unique F — v characteristic describing the dynamic 
behavior of the system is the S'-shape-like curve, and, 
accordingly, up- and down-switchings occur at the insta- 
bility points where \dv/dF\ =00. To derive rigorously 
this kind of the v(F) dependence a careful analysis in 
the vicinity of the critical point v c accounting for the in- 
teraction between the 2D and 3D modes and nonlinear 
pinning effects is needed. Leaving the detailed derivation 
for a forthcoming publication we point out that the re- 
lated scenario where the critical points merge and v(F) 
retains only the inflection point near v = v c is also pos- 
sible. In this degenerate case the plastic-elastic dynamic 
transition becomes non-hysteretic and can be detected by 
the position of the inflection point in v(F) dependence. 

Turning now to the case of weak disorder, V < /i 2 ; one 
can be easily verify, that in this case there will be always 
a non-zero solution for [i (note that for weak disorder 
Fj? D = V 4 /(i < F 2D — V 2 ). Combining this observation 
with the Eq. ( |Io| ) we conjecture the v — V diagram shown 
in Fig 2. 

Our discussion was restricted to T — 0, but the form 
of the F-v curve (S-shape, for example) is determinbcd 
by the intrinsic dynamic properties of the system and is 
to be stable with respect to thermal effects (as long as 
the the latter leave the periodic structure intact). An- 
other point is that although our consideration was fo- 
cused on an anisotropic CDW model, the form of the 
coarse-grained random force we used is not specific to 
CDWs but seem to be generic for any periodic structure 
subject to quenched disorder |L6 18|. We expect therefore 
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that the above ideas apply to a general case of periodic 
media driven through quenched disorder at finite tem- 
peratures. In particular, one can view the "strong depin- 
ning" switching scenario as the "zero-temperature pro- 
jection" of the sequential depinning or non-equilibrium 
freezing transition in the vortex system subject to strong 
disorder proposed in || . The instability force transcribes 
into a freezing force introduced for a moving vortex lat- 
tice and the independent layered motion of 2D CDWs 
maps onto a regime of plastic flow at the intermediate 
currents j c < j < jf (see inset in the Fig la). The 
details of the finite temperature behavior as well as the 
quantitative transcription of the ideas developed for the 
anisotropic CDW onto a general case of periodic driven 
media will be presented elsewhere. 

In conclusion, we described the depinning behavior in 
the driven periodic structures using the model of the lay- 
ered CDW as a prototype system. We have found that in 
case of strong disorder where the static state is decoupled 
the depinning occurs via a two stage process. First the 
pinned system experiences a continuous depinning into a 
plastically sliding state and at higher drives undergoes a 
second sharp hysteretic transition into a coherently mov- 
ing 3D state. The second transition is identified with 
the freezing transition proposed in S . For weak disorder 
depinning occurs into a coupled state. 
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FIG. 1. v — F transport characteristic for the strong pin- 
ning at T = 0. The dashed line displays viscous behavior 
in the absence of pinning. Inset: Dynamic phase diagram 
for periodic medium driven through strong disorder. The 
dashed line denotes the critical depinning force jc(T) at which 
crossover from creep plastic flow takes place. The solid line 
denotes the switching-like freezing transition at jf(T). 
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FIG. 2. v — V phase diagram: the system recovers coupling 
at v > vJV), the large velocity dependence of v c (V) is given 
by Eq.(px|)o 
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